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O ; Abstract 

We study the statistical properties of recurrence times in the self-excited Hawkes conditional 
Poisson process, the simplest extension of the Poisson process that takes into account how the past 
events influence the occurrence of future events. Specifically, we analyze the impact of the power 
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law distribution of fertilities with exponent a, where the fertility of an event is the number of 
aftershocks of first generation that it triggers, on the probability distribution function (pdf) f(r) 
of the recurrence times r between successive events. The other input of the model is an exponential 
Omori law quantifying the pdf of waiting times between an event and its first generation aftershocks, 
whose characteristic time scale is taken as our time unit. At short time scales, we discover two 
intermediate power law asymptotics, f(r) ~ 7- - ( 2 ~ a ) for r ^ r c and f{r) ~ T~ a for t c < t < 1, 
where r c is associated with the self-excited cascades of aftershocks. For 1 < t < l/i/, we find a 
constant plateau f{r) ~ const, while at long times, 1/v < r, /(t) ~ e~ UT has an exponential tail 
controlled by the arrival rate v of exogenous events. These results demonstrate a novel mechanism 
for the generation of power laws in the distribution of recurrence times, which results from a power 
law distribution of fertilities in the presence of self-excitation and cascades of triggering. 
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I. INTRODUCTION 



The statistics of 'recurrence times', defined as the random durations of time intervals 
between two consequent events, is widely used to characterize systems punctuated by short- 
duration occurrences interspersed between quiet phases. For instance, the statistics of re- 
currence times between earthquakes is the basis for hazard assessment in seismology The 
statistics of recurrence times has recently been the focus of researchers interested in the 



properties of different natural 
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25] and social systems 



The study of recurrence between earthquakes is perhaps the most advanced quantitatively 
due to the availability of data and the high involved stakes. The statistics of earthquake 
recurrence times in large geographic domains have been reported to be characterized by 
universal intermediate power law asymptotics, both for single homogeneous regions (3), 7| 
and when averaged over multiple regions 3|, |8| . These intermediate power laws, as well as the 
scaling properties of the distribution of recurrence times, were theoretically explained by the 
present authors 20. |21|. using the parsimonious ETAS model of earthquakes triggering 18], 
which is presently the benchmark model in statistical seismology. We recall that the acronym 
ETAS stands for Epidemic- Type Aftershocks and the ETAS model is an incarnation of the 
Hawkes self-excited conditional Poisson process 12Nl5|. Our previous works 2C| |2l| have 
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shown that one does not need to invoke new universal laws or fancy scaling in order to explain 
quantitatively with high accuracy the previously reported scaling laws , IsJ . In other 

words, the finding reported in jsljQ,^] do not contain evidence for any new physics/geophysics 
laws but constitute just a reformulation of the following laws: 



Earthquakes tend to trigger other earthquakes according to the same triggering mech- 
anism, independently of their magnitudes. 

The Omori-Utsu law for aftershocks, generalized into the phenomenon of earthquake 
triggering where all earthquakes are treated on the same footing, states that the rate 
of events that are triggered by a preceding event that occurred at time decays as 



hit) 



o<e<. 1, 



to > 0, 



t > 



(1) 



(to + t) 1+e ' 

The function fi(t) can also be interpreted as the probability density function (pdf) 
of the durations of the waiting time intervals between the reference "mother" event 
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and the triggered events of first generation corresponding to direct triggered by the 
mother event. The constant to describes a characteristic microscopic time scale of the 
generalized Omori law that ensures regularization at small time and normalization. 

Our previous analytical derivations 
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2l|, found in excellent agreement with empirical 



was essentially based on the long- memory of the Omori law (JTJ) , fi(t) ~ t^ 1 ^ 9 
for large t with < 9 < 1. However, it did not take into account the impact of heterogenous 
fertilities, which come in wildly varying values. Indeed, the number of daughters triggered by 
an earthquake of a given magnitude grows exponentially with its magnitude. For instance, a 
magnitude 8-earthquake may have tens of thousands of aftershocks of magnitude larger than 
2 while a magnitude 2-earthquake may generate no more than 0.1 earthquake on average 
of magnitude larger than 2 [16j]. Given the fact that the distribution of magnitudes is itself 
an exponentially decaying function of magnitudes (called the Gutenberg- Richter law), this 
translates into a heavy tail distribution of fertilities jl9], i.e. the distribution of the number 
of first generation events triggered by a given event has the following power law asymptotic: 

Pl (r) ~ r~ a -\ r ^ oo, a G (1, 2) . (2) 

Precisely, pi(r) is the probability that the random number Ri of first generation aftershocks 
triggered independently by a given mother event is equal to a given integer r. 



In fact, the main approximation in our previous work [20] , l2lj was to consider that, for 
the estimation of the distribution of recurrence times, it is sufficient to assume that each 
mother event triggers at most one event, so that the power law ([2]) is completely irrelevant. 
This surprising approximation was justified by the focus on the tail of the distribution of 
recurrence times, for which typically only one event, among the set of events triggered by a 
given earthquake, does contribute. 

The goal of the present paper is to reexamine this approximation and present an exact 
analysis of the impact of the power law form ([2]) of the distribution of fertilities on the 
distribution of the recurrence times. To make the analysis feasible and exact, we consider 
the case where the Omori law is no more heavy-tailed but has a shorter memory in the form 
of an exponential distribution, expressed by a suitable choice of time units in the form 

h(t) = e-*. (3) 

In addition to getting exact analytical expressions, our study of the case of an exponential 
memory kernel ([3]) is motivated by the fact that, for many applications, this is the default 



assumption 



0, y, y, 



10l . I 111 . |23j. Therefore, this parameterization has an genuine interest 
and intrinsic value. This exponential memory function should not be confused with the 
Poisson model, which has no memory. In contrast, the Hawkes process takes into account 
the full set of interactions between all past events and the future events, mediated by the 
influence function given by f\(t). 

The main result of the present paper is the exact expression for the full distribution f(t) 
of recurrence times that results from all the possible cascades of triggering of events over 
all generations. We make explicit the substantial dependence of the power law exponents 
characterizing the distribution f(t) on the exponent a of the power law tail (T5]) for the 
distribution of fertilities and on the branching ratio n, defined as the mean number of events 
of first generation triggered per event: 

n = E[R 1 ] . (4) 

The paper is organized as follows. Section 2 present the self-excited conditional Hawkes 
Poisson process and the main exact equations obtained using generating probability func- 
tions. Section 3 studies the probability of quiescence within a given fixed time interval 
and derives different statistical properties of earthquake clusters, such as the mean number 
of earthquakes in clusters and their mean duration. Section 4 presents the main results 
concerning the probability density functions (pdf ) of the waiting times between successive 
earthquakes. Section 5 summarizes our main results and concludes. The Appendix makes 
more specific the analytical model used in our derivations and describes useful statistical 
properties of first generation aftershocks. While we use the language of seismology and 
events are named 'earthquakes', our results obviously apply to the many natural and social- 
economic-financial systems in which self-excitation occurs. 



II. BRANCHING MODEL OF EARTHQUAKE TRIGGERING 

Before discussing the statistical properties of recurrence times, let us develop the statis- 
tical description of the random number of earthquakes occurring within the time window 
(t,t + r). The Hawkes model that we consider assumes that there are exogenous events 
(called "immigrants" in the literature on branching processes, or "noise events"), occurring 
spontaneously according to a Poissonian stationary flow statistics. Thus, successive instants 
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• • • < t-i < to < t\ < ti < . . . of the noise earthquakes belong to the stationary Poisson 
point process with mean rate v = const. Then, the Hawkes model assumes that any given 
noise earthquake, occurring at {tk}, triggers a total of R\ first generation earthquakes. We 
assume that the set of the total numbers {Ri} of first generation aftershocks triggered by 
noise earthquakes are iid random integers, with the same generating probability function 
(GPF) G\{z). Furthermore, the distribution of waiting times t between the time tk of a 
given noise earthquake and the occurrences of its aftershocks is assumed to be fx(t) defined 
by expression (|3]). In turn, each of the aftershocks of first generation triggered by a given 
noise earthquake also triggers independently its own first generation aftershocks, with the 
same statistical properties (same GPF) as the first generation aftershocks. 

Specifically, the Hawkes self-excited conditional Poisson process is defined by the following 
form of the Intensity function 

X(t\H t , V) = A noise (t) + Yl /(* - *<) ' ( 5 ) 

i\ti<t 

where the history H t = {tj}i<j<i tj t it <t<t it+1 includes all events that occurred before the 
present time t and the sum in expression (jSJ) runs over all past triggered events. The set 
of parameters is denoted by the symbol V. The term A no i se (t) means that there are some 
external noise (or immigrant, exogenous) sources occurring according to a Poisson process 
with intensity A no i se (i), which may be a function of time, but all other events can be both 
triggered by previous events and can themselves trigger their offsprings. This gives rise to 
the existence of many generations of events. In the sequel, we will consider only the case 
where A no i se (t) = A no i se is constant. 

Under the above definitions and assumptions, one can prove that the GPF of the total 
number of all earthquakes (including noise and triggered events of all generations) occurring 
within the time window (t, t + r) can be expressed as a product of two terms: 

e(z;r) = e out (z;r)-e in (z;r) . (6) 

The first term O ou t{z; r) is the GPF of the number of all earthquakes in (t, t + r) triggered 
by noise and triggered earthquakes that occurred up to time t. The second term Qi n (z; r) 
is the GPF of the number of all noise earthquakes that occurred within (t, t + r) and of all 
earthquakes triggered in that window by events also in (t, t+r). The factorization of Q(z; r) 
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given by expression ([6]) simply expresses the independence between the branching processes 
starting outside and within the time window (t, t + r). 



We have previously shown [2l|, |22j) that O ou t(- 2 ; T ) an d <dm( z ', T ) are given respectively by 
ou tO;r) =exp^/^ [0(*;*,r) - l]cftj , (7) 

and 

in (^r)=exp^ [zG(z;r)-l]d^ , (8) 
where the functions t, r) and G(z; r) satisfy the following nonlinear integral equations: 

G(z;t)=Q[H(z;t)) , 
H(z;T) = p(t)-zf 1 (r)®G(z;T) , (9) 
G(z;0) = l, #(>;0) = 0. 

and 

0(z;t,r) = Q[ft(z;t,r)] , 
H(z; t, t) = p(t + t) - g(z; t, r) ® f x {t) - zG(z; r) ® h{t + r) , (10) 
G(z; 0, r) = G(z; r), 0, r) = r) . 

The symbol £g> represents the convolution operator with respect to the repeating time argu- 
ments t or t. We have introduced the auxiliary function 

Q(y) := GiO- - y) , (11) 
and the cumulative distribution function (cdf) of the first generation aftershocks instants 

P(t)= tfi(t')dt'. (12) 



The functions G(z;t) and G(z;t,r) have an intuitive geometric sense, illustrated by 
figure [TJ G(z\ r) is the GPF of the random number of aftershocks of all generations triggered 
till the current time t = t by some noise earthquake (or some aftershock) that occurred at the 
origin of time t = 0. In turn, Q(z; t, r) is the GPF of the random number of all generations 
aftershocks triggered within the window (t, t + r) (for t > 0) by some noise earthquake that 
occurred at the origin of time. The GPFs G(z; r) and Q{z\ t, r) are related by their definition 
as follows: 

G(z;0,t) = G(z;t). (13) 
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FIG. 1: Schematic illustration of the geometric sense of function Q(z;t,r). It is the GPF of the 
random number R(t, r) of all generations aftershocks within the window (t, t + r) triggered by a 
noise earthquake that occurred at the origin of time t = (large solid arrow). The aftershocks 
within (t, t + r) are depicted by the small solid arrows. The aftershocks triggered outside the 
window (t, t + r) are depicted by the small dotted arrows. In the picture, R(t, r) = 5. 

A key input of the theory within the formalism of GPF is thus G\(y) or Q(y) (via its 
definition (fTTjl ). Given expression (j2J), we show in the Appendix that a convenient form is 

Gx(z) = l- n(l- z) + k(1- z) a , ae(l,2], (14) 

so that the corresponding auxiliary function Q(y) ( fill takes the form 

Q(y) = l-ny + Ky a . (15) 

III. PROBABILITY OF QUIESCENCE, MEAN NUMBER OF EARTHQUAKES 
AND MEAN DURATION OF CLUSTERS 

A. Probability of quiescence 

For z = 0, relation (jHJ) reduces to 

P(r) = Pout(r)-P in (r), (16) 
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where 

P(r) = &(z = 0;r) (17) 

is the probability that there are no earthquakes (including the noise earthquakes and their 
aftershocks of all generations) within the window (t, t + r). P(t) can be decomposed as the 
product of the probability Pj n (r) that no noise earthquakes occur within (t,t + r) and the 
probability P ut(' 7 ") that no aftershocks occur within (t,t + r) that could have been triggered 
by noise earthquakes and their aftershocks that occurred before and until time t. 
In the general case, Pi n (r) is given by 

Pin(r) = e~ vr . (18) 

In turn, P ou t( r ) = Qout( z = 0; r) by definition. Using relation ([7]) and equations Q), ffTUl) . 
we obtain 

P out (r) = exp(i/jf [g(t,r)-l]dtj , (19) 



where 



is defined by 



g(t,r) = g(z = 0;t,r) (20) 



g(t,r) = Q[H(t,r)} . (21) 
The auxiliary function 7i(t, r) is solution of the nonlinear integral equation 

H(t,T)= P {t + T)-Q[H{t,T)]®f 1 {t) , 

(22) 

U(0,t)= P (t) . 

B. Solution of equation ([22]) and determination of P(r) for the exponential pdf 

hit) 

Using the exponential form of the Omori law ([3]), it is possible to obtain an exact analytical 
solution of equation (j22l) . which gives us the possibility to explore in detail the probabilistic 
properties of recurrence times. 

Using the form ([3]), it is easy to show that equation ( [22]) reduces to the initial value 
problem 

— + n + Q[H} = l, H(0,t) = p(t). (23) 
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Using expression (jT5l) for Q(y) leads to 



— + (1 - n)H + k7T = 0, H(Q,t)= P (t) 



whose solution is given by 



H(t,r)= [(1-e 



where 



-r-jl-a e (a-l)(l-»)t + ^ / g (o-l)(l-»»)t _ jjjVM) 
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1 — n 



We have used the fact that 

p(r) = 1 _ e -r , 

as derived from the exponential pdf fi(t) given by ([3]) and definition (|12p . 
We can now rewrite probability P ou t(T) f jT9T) in the form 



PoutW = e 



-i/F(r) 



where 



[l-^(t,r)]dt 



Taking into account expression f[25|) and the equality 



0(t, r) = Q[H(t, T)]=l-nH + kH c 
the explicit calculation of integral fl29|) yields 



F(r) = F(n, /c,a,p) 
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V(i-") 



7P 



1 a - 2 



(a — 1)(1 — n)|_ \p + 7p a ' a — 1 ' a — 1 



1 -nVB 



7P C 



a 



p + 7p Q ' a — 1 ' 1 — a 
where p = p(r) and B(x; a, b) is the incomplete beta function 



B(x; a, b) 



1 - s) h - l ds . 



(24) 

(25) 

(26) 
(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



In view of the key role played by the function F(n, k, a, p) in the following, it is useful to 
describe some of its properties. A first result of interest is its limit behavior as the branching 
ratio n tends to 1. Recall that this limit corresponds to the critical regime of the Hawkes 



process, separating the subcritical phase n < 1 and the supercritical phase n > 1. For 
n < 1, each noise earthquake has only a finite number of aftershocks. For n > 1, there is 
a non-zero probability that a single noise earthquake may generate an infinite number and 
infinitely long-lives sequence of aftershocks. The relevant physical regime is thus n < 1 and 
the boundary value n = 1 plays a special role, especially when one remembers that n can 
also be interpreted as the ratio of the total number of triggered events to the total number of 
events [17J]. Hence, when n — >■ 1, most of the observed activity is endogenous, i.e., triggered 
by past activity. 

Thus, the limit of F(n, k, a, p) as n — > 1 reads 



2-a 



F(k, a, p) = lim F(n, k, a, p) = 

n-tl K,{2 — a) 



- p, 1 < a < 2. 



(33) 



For n < 1, it is convenient to choose values of a that take the form 

1 



a = 1 + 



ni 



(34) 



so that it is possible to express the function 



J* m (n, k, p) = F in, k, H , p 

m 



(35) 



under the form 



rap 



1 — n 



m , 



m—1 



In 1 



l/m 



ln r , 



-P 



l/m 



(36) 



1 — 71 / \ 1—71 

The auxiliary function ln m (l + x) is defined as the sum of the first m terms of the Taylor 
series expansion of the logarithm function m(l + x) with respect to x: 

\k 



ln m (l + x) = - ^ 



k=l 



k 



(37) 



For m = 1 (a = 2), we have 



1 + 



np 
1 — n 



P ■ 



JVn, k, p) — — In 

K 

For m = 2 (a = 3/2), we have 

2 2 f k fo\ 

T 2 {n,K,p) = - y /p-p-—{l- n ) In (^1 + ^3^) 



(38) 



(39) 
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For m = 3 (a = 4/3), we have 



K 2k 
-j(l -„)•]» 1 + ^- 



C. Mean duration of seismic clusters 



i. Seismic clusters of all types 

Let us consider the random duration A& of the aftershock cluster triggered by the fcth 
noise earthquake that occurred at time t^. By definition, 

Afc = tLt ~ t k , (41) 

where tf ast is the occurrence time of the last of its triggered aftershocks over all generations. 
The mean value of A^ is by definition (A) = J °° g ■ w(g)dg, where w(g) is the pdf of the 
random cluster durations {A k }. 

By hypothesis, the instants {tk} of the noise earthquakes form a Poissonian point process 
with mean rate v. Moreover, within the Hawkes branching process model, the clusters 
durations {A^} are iid random variables. One can easily show that the probability P ou t(T) 
given by ( 1281) . that no aftershocks occur within (t, t + r) that could have been triggered 
by noise earthquakes and their aftershocks that occurred before and until time t, take the 
following value in the limit r — > +oo: 

P out (oo) = e-'< A > , (42) 

Thus, P out (oo) is the probability that all noise earthquakes that occurred up to time t do 
not trigger any aftershock after t. It follows from relation ff28l) and fj4"2l that 

(A) = Foo(n, K,a, p = 1) . (43) 

Now, F oo(n, n,a, p = 1) can be obtained as the value of the function F(t) given by fl3T|) 
at r — > +oo: 

Foo(n, k, a, p) = F(oo) . (44) 
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FIG. 2: Dependence with respect to the exponent a of the probability P ou t(oo) (|42|) that no 
aftershocks occur within (i, t + r) that could have been triggered by noise earthquakes and their 
aftershocks that occurred before and until time t. The plot corresponds to the critical case n = 1 
with k = 0.25. Bottom to top: v = 0.1; 0.05; 0.01; 0.001. 



We obtain 



F(n, K,a, p = 1) 



7 



l/(l-a) 



(a - l)(l-n) 
fl - n)5 



K 



1 a-2 



k + 1 — n'a — l'a — 1 



1 



^ + 1 — n a — 1 1 — a 
In particular, in the critical case n = 1, we have 

1 



F(l, k, a, V 



k(2 -a) 



and thus 



(A) 



1 



- 1, n = 1, 



(45) 



(46) 



(47) 



k(2 - a) 

Figure [2] shows the probability P out (oo) ( )42l) as a function of the exponent a, in the 
critical case n — 1. 
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2. Seismic clusters with at least m > 1 aftershock 

The mean duration (A) of clusters given by expression ( 143|) with f H5|) includes the contri- 
bution of the empty clusters for which the noise earthquake does not trigger any aftershock. 
It is thus interesting to evaluate another derived quantity (A 1 ), defined as the mean du- 
ration of clusters that contain at least one aftershock. To get (A 1 ), we divide (A) by the 
probability that the number of aftershocks is strictly positive, 

(A 1 ) = (A)/Pr{i? > 0} , (48) 

where Py{R > 0} is the probability that the number R of aftershocks is positive. Accord- 
ingly, one introduce the mean rate v l of the non-empty clusters equal to 



v 



u/Pr{R>0}, (49) 



where v is the mean rate of noise earthquakes. 
Obviously, Pi{R > 0} is given by 

Pr{R > 0} = 1 -pi(0) , (50) 

where pi(0) is the probability that a noise earthquake does not trigger any first generation 
aftershock at all. Using the parameterization defined in the Appendix for the Hawkes model, 
Pi(0) is given by expression ( 185]) . leading to 



v 1 = v ■ (n — k) 



and (A 1 ) = F - (n '"' a ' 1} . (51) 



n — k 



In particular, in the critical case n — 1, the mean duration of the non-empty clusters is given 
by 

x l-/c(2-q) 
(A } = (2-a)«(l-«) ■ (52) 

As an example, taking k = 0.25 and a = 1.5 yields a mean duration of the non-empty clusters 
in the critical case equal to (A 1 ) ~ 9.33. Recall that the unit time is the characteristic decay 
time of the Omori law fi(t) (J3j). 

This allows us to define regimes of low seismicity as characterized by the following in- 
equalities 

v ■ (A) < 1 <=> v 1 - (A 1 ) < 1 , (53) 
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which means that clusters are well individualized, being separated by comparatively long 
quiet time intervals. 

It is useful to generalize the mean duration (A 1 ) of clusters that contain at least one 
aftershock to the mean durations (A m ) of the clusters that contain at least m aftershocks. 
We now derive the general equation allowing one to calculate these (A m )'s. Using the total 
probability formula, one can represent the pdf w(g) of clusters durations in the form 

oo 

w{g) = p(0)6(g) + J2p(j)w(g\j) , (54) 
j'=i 

where p(j) is the probability that a given noise earthquake triggers j aftershocks of all 
generations, and w(g\j) is the conditional pdf of cluster durations under the condition that 
the number of aftershocks is equal to j. Accordingly, the pdf w(g\j ^ m) of the durations 
of the clusters that have m or more aftershocks is equal to 

m—1 

w(g) - ^2p(])w(g\j) 
w(g\j > m) = 3 -^ I{ . (55) 

j=Q 

The corresponding conditional expectation (A m ) is equal to 

m—1 



(A)- J>0')(A|j> 



(A m ) = / g w(g\j > m)dg = ^- . (56) 



3=0 

In particular, taking into account that 

p(0) = Pi(0) = 1 — n + k, p(l) = Pi(l) = n — an, 
w(g\l) = fi(g) =}► (A|j) = l, 



(57) 



we obtain 



(A 2 ) - (^-Pt 1 ) m 

{A) -i- P (o)-p(iy (58) 

Using (|57|) for n = 1 and the expression (|47|) for (A), we obtain the mean duration of clusters 
containing more than one aftershock: 

2 I - k(2 - a)(2 - an) 

(A } = «?(a-l)(2-a) ■ (59) 
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FIG. 3: Dependences of the cluster durations (A) (|47p (bottom curve), (A 1 ) ([52]) (middle curve) 
and (A 2 ) (5J) (top curve) as functions of the exponent a. The plot corresponds to the critical case 
n = 1 with re = 0.25. Recall that (A m ) is defined as the mean duration of the clusters that contain 
at least m aftershocks. In our notations, (A) corresponds formally to (A ), i.e., it also takes into 
account the empty clusters for which the noise earthquake does not trigger any aftershock. 

The dependences of the cluster durations (A) (HTj) . (A 1 ) (]5"2]) and (A 2 ) (]5"§]) as functions of 
the exponent a are shown in figure Note the large jump in mean durations of clusters 
containing at least two aftershocks compared with clusters containing at least one aftershock. 



IV. PDF OF RECURRENCE INTERVALS 



A. General relations 

The knowledge of the exact probability P(r) f[T6]) allows one to calculate exactly the pdf 
f(r) of the random waiting times {T k } between subsequent earthquakes. Indeed, a general 
result of the theory of point processes states that 

d 2 P(r) 

m = ^y-j^ 1 > (6°) 
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where (r) = E [Tf-\ denotes the mean waiting time between subsequent earthquakes. There- 
fore, the complementary cumulative distribution function (ccdf ) of the random waiting times 
is equal to 

*(r) = Pr {T > r} = -(r)^^ . (61) 

dr 



By normalization, ^(0) = 1, so that 

1 dP(r) 



(t) dr 

Using expressions ffTB]) . f fT8"j) and (}2"8"]) . we have 



(62) 



T=0 



P(r) = e -^W-^ . (63) 



Making explicit F(r) with expression (I3"T|) yields 

dP(r) = ^ l-np(r)+Kp"(r) ^ (T 
dr 1 - n + Kp a ~ l {T) 

Using ( )62|) . we have 



F(t)-ut 



(64) 



. . 1 — n T . . 1 — n dPir) 

r= => tf(r = 65 

v v dr 



and finally obtain 



*(r) = (1 - n) l - np{T) + Kp ^ e-^- . (66) 

Differentiating this last expression ( )66|) with respect to r yields the pdf /(r) of waiting times 
between successive earthquakes: 

/(r) = ^{n,K,a,p{r),v) , (67) 



where 



and 



$(n, k, a, p, z/) = 
[„4(n, «, a, p) + ^ • B(n, k, a, p)\ e --{n^ a , p)+T ) ^ 

A(n, k, a, p) = (1 — n)(l — p) x 
- n) + «(a - 1 + (2n - a)p)p Q ^ 2 - K 2 p 2(a_1) 



(68) 



1 - n + Kp a - X ) 2 

2 

^ I I — ; ) n -+- I-: li ' \ 

B(n, k, a, p) = (1 — n) 



(69) 



1 — np + np c 



1 — n + Kp 



,a-l 
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Recall that p represents pit), which is defined by expression (fT2|) . 

In the following subsections, we analyze in detail expressions fl68|) and (159]) in order to 
derive the properties of the pdf /(r) (1571) . Expression (1551) suggests that it is natural to 
decompose the analysis of /(r) into two discussions, one centered on the term surviving in 
the limit v — > and the other one. The two next subsections analyze these two terms in 
turn. 



B. Case v — > 

Taking the limit v — > amounts to neglecting the occurrence of any noise earthquake 
within the window (t, t + r) of analysis. As shown in the next subsection, this first case 
already reveals interesting properties of the pdf /(r), which remain valid in the general case 
v > 0. 

Putting v = in expression ( 1551) and using ( 1571) . we have 

/(r) = ^[n,/c,a,p(r)] (*> = 0) , (70) 

where the function A(n, k, a, p) is given by expression (15U1) . 

The main asymptotics of the function A(n, k, a, p) are respectively 

• at p <C 1: 

a-2 

^(n,K,a,p)~K(l-n)(a-l) H , p<l. (71) 

(1 — n + ftp J 

This regime p 1 corresponds to r < 1 and thus p ~ r. Relation ( ITT]) thus leads to 

/(r)^«(l-n)(a-l) ^ ] r«l. (72) 

(1 — n + kt L y 

• At p -> 1: 

^(n,«,a,p)~C(l-p) , 1-P«1, C= ( 1 ~ ra )^- K ) , (73) 

1 — n + k 

This second asymptotic p — >■ 1 corresponds to 1 — p <C 1 , which is equivalent to the 
condition r ^> 1. Using expression (1271) . relation (173|) leads to 

/(r)~Ce- r , r>l. (74) 
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The denominator of expression (!72|) determines a new characteristic time scale 

(75) 



i 

1 — 



K 

from which one can define a critical value for the branching ratio 



1 - k , (76) 



such that t c > 1 for n < n c and r c < 1 for n > n c . We shall refer to the case n < n c as the 
subcritical regime, while n c < n < 1 is called the near-critical regime. 

For n < n c (subcritical regime), r c > 1, and one may replace the asymptotics (l7"2l) by the 
pure power law: 

f(r) ~ k ^ , n < n c . (77) 



In contrast, in the near-critical case n c < n < 1, expression (J72J) leads to two power 
asymptotics: 

[ K r^ r ~ (2 ~ a) > r<<Tc ' 

f(r) ~{" n c <n<\. (78) 

a -I 

1-n r- Q , r c <r<l, 

v k 

Figure H] shows the pdf f(r) (1701) of the waiting times between successive earthquakes for 
v = 0, a = 1.5, k = 0.25 (n c = 0.75), and in the near-critical case n = 0.999. The two power 
law asymptotics (ITBl) and the exponential asymptotics fTM|) for large r's are clearly visible. 
Figure [5] is the same as figure HJ except for the value n = 0.9. Although, formally, this 
value belongs also to the near-critical case (n = 0.9 > n c = 0.75), the intermediate power 
law asymptotics r~ a is barely visible and, for any r C 1, the subcritical power asymptotics 
r -( 2 -") dominates at short times. 

Figure [6] shows log-log plots of the pdf /(r) for v = 0, k = 0.25, a = 1.5 and different 
values of the branching ratio n belonging to the subcritical regime. The subcritical power 
law asymptotics r~( 2-a ) is dominant and the different pdf's are similar. 



C. General case v ^ 

We now take into account the occurrence of noise earthquakes within the window (t,t + r) 
of analysis. The main qualitative difference between this case and the previous one v = is 
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FIG. 4: Solid line: Plot of the pdf f(r) (|70p of waiting times between successive earthquakes for 
v = 0, k = 0.25, a = 1.5 and n = 0.999. The dotted lines show the power law asymptotics ([751) 
and the exponential asymptotic behavior (|74p . 




FIG. 5: Same as figure H] except for the value of the branching ratio n = 0.9. 
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10" 6 10" 4 10" 2 10° 

r 

FIG. 6: Solid line: Log-log plot of the pdf /(t) (|70p for v = 0, n = 0.25, a = 1.5 and for three n 
values in the subcritical regime. Top to bottom: n = 0.8; 0.6; 0.4. The dotted straight line is the 
subcritical power law asymptotics (f77|) . 

the replacement of the large r asymptotics ( 1741) by 

/(r) ~ (1 - n) ( n ~ K e -r + v \ e -K<A>+r) j r » j . (79) 

The interesting regime of well-defined clusters occurs for v <C 1, for which the typical 
waiting time between noise earthquakes is much larger than the characteristic decay time of 
the Omori law, which is the typical waiting time between a noise earthquake and its after- 
shocks. In the interval 1<t<1/i/, expression (I7§|) reduces to (1741 and this exponential 
decay can be replaced by the plateau 

/(r) ~ (1 - n)u , Kr<l/v. (80) 

For r > 1/V, expression (I7§1) simplifies into 

/(r) ~ (1 - n)z/e-^ A > e— , r > l/v . (81) 

These different regimes are illustrated in figure [7J which depicts the pdf of waiting times 
between successive earthquakes for k = 0.25, a = 1.5, v = 0.999 and three values of the 
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FIG. 7: Solid line: Loglog plots of the recurrence times pdf /(r) ([67]) in the case n = 0.25, a = 1.5 
n = 0.999 and three subcritical v values: v = 0.1; 10~ 3 and v = 10~ 5 . Dotted straight lines are the 
subcritical and near-critical power asymptotics (J73|) 

positive parameter z/ = 10 _1 ; 10~ 3 ; 10~ 5 . One can clearly observe the subcritical and near- 
critical power law asymptotics ( l78l) . as well as the plateau ( )80|) joining the exponential 
asymptotics (TTijl and ( 18TT) . The plateau is especially visible for the smallest values v = 10~ 3 
and v = 10~ 5 . 



V. CONCLUSION 



In this paper, we have studied the impact of the power law form (T5J) of the distribution 
of fertilities on the distribution of the recurrence times. In contrast with previous studies, 
we have considered the simple case where the Omori law is not heavy-tailed but is described 
by an exponential function. Motivated by real applications, this choice allows us to develop 
an exact analytical treatment using the formalism of generating probability functions. 

Our analysis emphasizes the importance of three time scales controlling the different 
regimes of the probability density function (pdf) /(r) of waiting times between successive 
earthquakes: 
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• the characteristic decay time of the exponential Omori law fi(t) ([3]) describing the 
occurrence of first generation aftershocks, taken as our time unit, 

• the average waiting time 1/v between two successive "noise earthquakes" , which con- 
stitute the exogenous sources of the self-excited processes followed by their aftershocks, 

• the characteristic time r c defined by expression (|75p associated with the self-excited 
cascades of aftershocks. 

In the interesting and relevant regime where earthquake clusters are well defined, namely 
when the typical waiting time till the first aftershocks (unit time) is smaller than the waiting 
times 1/v between noise earthquakes, we have found that the pdf of recurrence times exhibits 
several intermediate power law asymptotics: 

1. For t < t c , f(r) ~ T -( 2 - a \ 

2. For r c < r < 1, /(r) ~ T~ a . 

3. For l«r«l/i/, /(r) ~ (1 — n)v = const. 

4. For l/v<r, f(r) ~ e~ u ^ e~ UT . 

In these formulas, a is the exponent of the power law distribution of fertilities px(r) ~ r ~ a -^ 
which is the pdf of the number of first generation aftershocks triggered by a given event 
of any type. In turn, n stands for the branching ratio defined by equation (j3J), i.e. the 
average number of daughters of first generation per mother event, and (A) is the mean 
duration of a cluster that starts with a noise earthquake and ends with its last aftershock 
over all generation. It is given by expression (j43j) . 

Only the first two intermediate asymptotics /(r) ~ t~( 2-q ) and /(r) ~ r~ a at short time 
scales reflect the influence of the power law distribution of fertilities (j2J), which is revealed 
by the remarkable effect of the cascade of triggering over the population of aftershocks of 
many different generations. 



Final 
work 20 



y, let us stress the differences between the present investigation and our previous 



2l| on the same problem. In Refs. 20|, |2l| , we determined the asymptotic behavior 



at long times of the distribution of recurrence times, under the approximation that it was 
sufficient to consider only one aftershock at most per mother event ('noise earthquake' in the 
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present terminology). In addition, we considered the standard power law Omori law ([T]) and 
not the exponential law ([3]). By performing a detailed analysis made possible by the use of 
the exactly tractable exponential Omori law ([3]), the present paper has thus demonstrated 
the existence of additional short-time intermediate asymptotics that reveal the distribution 
of fertilities. This opens the possibility to estimate the exponent a of the distribution of 
cluster sizes from purely dynamic measures of activity. 
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Appendix: Statistics of first generation aftershocks 



Given the power law ([2]) for the right tail of the pdf {pi(r)} of the number R\ of the 
first generation aftershocks probabilities, we show that the leading relevant terms of the 
expansion of the GPF Gi(z) in powers of (1 — z) take the form 

Gi(z) = l-n(l-z) + K(l-z) a , ae(l,2], (82) 

so that the corresponding auxiliary function Q(y) (TTTT) takes the form 

Q(y) = l-ny + Ky a . (83) 

They depend on the branching ratio n defined in (j3J), the exponent a of the power law 
distribution fl2]) of the number of first-generation aftershocks. The additional scale parameter 
k satisfies the following inequalities 

< oik < n, n ^ 1 , 

(84) 

n — 1 < aK < n, n > 1 , 
which ensure the necessary constraints 

0^Pi(0)<l, and 0<r(1)<1. 

Rather than deriving the form (1821) from (j2j), it is more convenient to show that the tail 
of the pdf {pi{r)} whose GPF is given by (1821) is the power law (121 . Given (I82p and the 
definition linking G\{z) to pi(r), namely G\(z) = Ylj=o Pi{ r ) zV ■> we obtain 

Pi(0) = 1 — n + k, — n — aK , (85) 

and 



pi(r) = k(-1) 



a\ K{-l) r r(a + i) 



(86) 



rj r(r + l)r(a - r + 1) 

r^2, a 6 (1,2). 
Using the properties of gamma functions and in particular the well-known equality 

r(*)r(i -z) = ^— , (87) 

sm it z 



we obtain that 



;_l)T(a-r+l) = = ? - y aG(l,2). (88) 

T(r — a) sm[7r(a — 1)J 
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Accordingly, expression (|86|) takes the form 



/ x T(r — a) k .. 

Pl(r) = c 'fFTTT c '^Way < 89) 



using the asymptotic relation 



r(r + a) a _ b 



T(r + b) 

we finally recover the power law (J2]). 

The case where the exponent a = 2 in expression (182]) of requires a special mention. 

Indeed, this form describes the special situation in which each noise earthquake (and any 
aftershock as well) can trigger not more than two first generation aftershocks. Accordingly, 
there are, in general, only three nonzero probabilities 

Pl (0) = l-n + /c, pi(l)=n-2«, pi(2) = « (a = 2) . (90) 

This special situation arises due to the fact that the GPF G±(z) has been truncated beyond 
the quadratic order (1 — z) 2 . It will not be considered further in this paper. 
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